Classical Investigation of Frustrated and Dimerized Heisenberg Chains 
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We have considered the ID dimerized frustrated antiferromagnetic (ferromagnetic) Heisenberg 
model with arbitrary spin S. The exact classical magnetic phase diagram at zero temperature is 
determined using the LK cluster method. Cluster method results, show that the classical ground 
state phase diagram of the model is very rich including first and second-order phase transitions. In 
the absence of the dimerization, a second-order phase transition occurs between antiferromagnetic 
(ferromagnetic) and spiral phases at the critical frustration a c — ±0.25. In the vicinity of the 
critical points a c , the exact classical critical exponent of the spiral order parameter is found 1/2. In 
the case of dimerized chain (5^0), the spiral order shows stability and exists in some part of the 
ground state phase diagram. We have found two first-order critical lines in the ground state phase 
diagram. These critical lines separate the antiferromagnetic from spiral phase. 

PACS numbers: 75.10.Pq, 75.10.Hk 
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I. INTRODUCTION 

During the last decades several classical techniques 
such as the well-known Luttinger-Tisza method^, vertex 
model 2 - and so on have been introduced to solve classi- 
cal Hamiltonian exactly. The Luttinger-Tisza method 
is more effective in systems with bilinear interactions 
and vertex model usually applied for treating frustrated 
modcl^. 

In a very recent work^, T. Kaplan have used a kind 
of cluster method, hereafter simply LK method, which 
is based on a block of three spins to solve frustrated 
classical Heisenberg model in one dimension with added 
nearest neighbor biquadratic exchange interactions. He 
asserted that the LK method is not limited to one dimen- 
sion or to translationally invariant spin Hamiltonians^ 
and expanded his approach to determine the phase dia- 
gram of frustrated classical Heisenberg and XY models 
with added nearest neighbor biquadratic exchange inter- 
actions in d = 2 dimension 7 -. In order to check the valid- 
ity of Kaplan's phase diagram conjecture, we have inves- 
tigated his models form quantum point of view for spin-i 
with an accurate algorithm (Lanczos method), and our 
results, which will be presented elsewhere^, showed that 
LK method, albeit is a classical approach but has the 
capability to work for aspects of a quantum treatment. 

Actually, these are our stimulating reasons to take a 
quite well known frustrated and dimerized Heisenberg 
model and determine its classical ground state phase dia- 
gram exactly with strong, but not well known LK method 
which is able to solve problems rigorously^. Let us start 
with definition of the dimerized and frustrated Heisen- 
berg model as follow 

H= Ji^[l + (-l)"<s]s„-S„ +1 + J 2 ^S„-S„+ 2 , (1) 

n n 

where S n is the nth classical vector of the length S. 
A spin system is frustrated when the global order be- 



cause of the competition of different kinds of interaction 
is incompatible with the local order, so chain with both 
antiferromagnetic-antifcrromagnetic exchanges (Ji > 
0, J2 > 0) and ferromagnetic-antifcrromagnctic ex- 
changes (Ji < 0, J2 > 0), hereafter simply AF-AF and 
F-AF respectively, are frustrated. 

Quantum study of this model is well done for the spin- 
1/2 and spin-1 chains&i2»iiii2 1 i2iiiiiLi2. It is found that 
the quantum fluctuations play a prominent role at zero 
temperature in the ground state phase diagram of the 
models. This model shows a dimerization transition at 
S = 8 C . For S = 1/2, the transition point is always 
S c — 0, since the Lieb-Schultz-Mattis theorem implies 
either gapless excitations or two-fold degeneracy of the 
ground states at 5 = 0. In fact, on the undimerized line 
S = 0, for AF-AF case, there exists a critical frustration 
parameter a c = 0.241 1 17 ' 18 . For a < a c the system is a 
gapless Tomonaga-Luttingcr Liquid (TLL); that is, the 
dimerization transition at 8 = S c = is of second or- 
der. On the other hand, for a > a c , the ground state 
is doubly degenerate, exhibiting a spontaneous dimeriza- 
tion. This implies a first-order dimerization transition at 
5 = 5 C = 0. For S = 1, S = (for a small a) belongs to 
Haldane phase and does not exhibit a transition line. In- 
stead, dimerization transition between the Haldane phase 
and the dimerized phase occur o 19 ' 20 at a finite S c , which 
depends on the frustration a. Albeit the shape of the 
phase diagram is thus different, topology of the phase di- 
agram is a bit similar to that for S = 1/2. In fact, also 
for S = 1, there is a critical frustration a c ; the transition 
is second order with the critical point described by a TLL 
for a < a c , and first order for a > a c . 

But there is not a classical clear picture of different 
ground state phases of the mentioned model. Having a 
classical picture, from one hand help us to know that 
quantum fluctuations destroy which one of the classical 
orderings. On the other hand for arbitrary large spin 
model, the classical picture is the same with the quan- 
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FIG. 1: (Color online.) The schematic picture of dimerized 
spins chain. 



turn picture of the ground state phase diagram. In this 
work we focus on the ID frustrated and dimerized sys- 
tems with arbitrary spin S (see FIG.[T]). To find the exact 
classical ground state phase diagram of the model, the 
LK cluster method is used. In the absence of the dimcr- 
ization, by increasing the frustration a classical phase 
transition occurs at a c = +0.25 (—0.25) from the anti- 
ferromagnetic (ferromagnetic) phase into the spiral mag- 
netic phase. Our results show that the dimerization pa- 
rameter induces new magnetic phases including stripe- 
antiferromagnetic phase (or uud and duu phases). Exist- 
ing of these magnetic phases is independent of length of 
spins. 

The outline of the paper is as follows. In forthcoming 
section we will extensively explain the LK method with 
implementing it to our model and in the section III we 
will summarize our results. 



II. THE LK CLUSTER METHOD 

In order to implement LK method we follow exactly 
the procedure in Rcf. [2]. Without losing the generality 
and setting periodic boundary conditions, Eq.([T]) can be 
rewritten as: 



H c — /J /i c (Si_i, Si, Si + i^ 



(2) 



where the "cluster energy" involve three neighboring 
spins is 

MSi,s 2 ,s 3 ) = Aj(i_£) Sl .s2 + (i + 5)S2-s 3 )} 

+ J 2 (Si-S 3 ). (3) 

It is clear that 

H c > y^mm h c (S j-x,S j,S j+i). (4) 

i 

To minimize h c respect spins directions, we first consider 
coplanar spins, and label the angles 9, 9' made by the end 
spins with the central spin (see FIG. [2]) which in coplanar 
case we set <j) = 0, <f)' = 0. The cluster energy is given by 

h c (9,9') = 5 2 {(l T ^)cos0+(i±^)cos0' 

+ acos(0-0')}, (5) 





FIG. 2: Spins cluster configuration. 



where a = J2/ J\ ■ Minimizing h c respect 
following equation: 



1' gives the 



dhc 
89 
dh c 



2 

2 



(l-5)sin0 + 2a!sin(0-0') 
(1 + 5) sm6' -2asin(0-0') 



(G) 



Let's first deal with a case without dimerization, by set- 
ting 5 = in Eq. ([5]) we have 



dh c 


S 2 


89 
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8h c 


52 r 


89' ~ 


2 



sin ( 



2a sin(9 — 



sm9' -2asin(9 -9') 



(7) 



its solutions are 



(0,0) ) (7r,7r),(7r,0),(0,7r), 
(<9 ,-6>o), (Spiral- type) 



where 



cos 9o = — — 
4a 



. 1 

H> 5 - 



(8) 



The solutions (ir, it), (0, 0) are related to collinear anti- 
ferromagnetic and ferromagnetic states respectively. The 
antiferromagnetic (ferromagnetic) state will minimize the 
energy in the case of J x > (J\ < 0). Solutions 
(7r,0) -> (4-, tit) an d (0,7r) (t,t, I) ar e degenerate 
states and show spins propagate in the down-up-up and 
up-up-down respectively^. Spiral state (9 ,—9 ) with 
uniform rotation is also degenerate state. In following 
we present results of the antiferromagnetic case J > 0. 

By setting minimization conditions into the Eq. ([5]) we 
have the following energies: 



^antiferro 
^uudd 

^spiral 



h c (w, 7T) 
h c (0,n) 

h c (9o, - 



S 2 ( 



5 2 ( 



l + o), 
a), 

8a 



a). (9) 



By equating these energies in pairs we have found only 
one critical point, a c = 0.25. Because of the continuity in 
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the derivative , a second-order phase transition occurs 
when passing trough a = 0.25. The ground state is in the 
antifcrromagnetic phase in the region of the frustration 
a < 0.25 and in the spiral phase in region a > 0.25. In 
general, the antiferromagnetic phase is recognized by the 
non-zero value of the Neel order parameter defined as 

m« = l£(-i)»s«, (io) 

n 

and the spiral phase in the ground state phase diagram 
of the spin systems is characterized by the nonzero value 
of the spiral order parameter 

x = ^£is„xs n+1 |. ( n ) 

n 

Using Eq. (JHJ we have found the spiral order parameter 
as 

X = , a < a c , 

* = S 2 \/l-^f , a>a c . (12) 

In FIG O we have plotted the spiral order parameter 
as a function of the frustration parameter a for the non- 
dimerized model (6 = 0). As is clearly seen from this 
figure, there is no long-range spiral order in the region 
of frustration a < a c = 0.25. However, in the region 
a > a c = 0.25 the spins of the system show a profound 
spiral order which grows by increasing the frustration 
parameter a. 

It has been discovered that continuous phase transi- 
tions have many interesting properties. The phenomena 
associated with continuous phase transitions are called 
critical phenomena, due to their association with critical 
points. It turns out that continuous phase transitions can 
be characterized by parameters known as critical expo- 
nents. Critical exponents describe the behavior of phys- 
ical quantities near continuous phase transitions. It is 
believed, that they are universal, i.e. they do not depend 
on the details of the physical system. Our analytical re- 
sults show that the spiral order parameter \ approaches 
zero in a singular fashion as a approaches a c , vanishing 
asymptotically as 

X«(l--) 1/2 , (13) 
a 

which shows that the critical exponent for spiral order 
parameter is a simple fraction e = 1/2. 

Now, we back to our original problem 5^0, find- 
ing the exact ground state phase diagram of the classical 
frustrated and dimerized Heisenberg chains. One can 
immediately see the possibility of having two spirals, one 
on the even sites, the other on the odd sites, both with 
the same wave length, but with a phase difference as de- 
scribed in the following: 




a 



FIG. 3: (Color online). The spiral order parameter as a 
function of the frustration parameter a in the case of non- 
dimerized chain (S — 0). 

To use the cluster approach for the case of non-zero 
dimerization, 5^0, one must consider two types of 
cluster, one being (• • •), the other being (• • •). 
Let us label the spins in the first cluster (Si S2S3), 
and the second (52 5354) and continue (53 545s) , 
(54 S5S6), ■ • ■ • The 2-spiral form assumes a simple spi- 
ral on the even sites, and a simple spiral on odd sites, 
both with the same turn angle between spins, which we 
called 26*o. Calling the angles — 6q + e and #0 + e, respec- 
tively, of Si and 53 with the center spin 52 in the first 
cluster, then gives 2# as the angle between Si and S3. 
Then, preserving the angle between S2 and S3 in the next 
cluster, (which is now — 8q — e , since the central spin is 
now S3) and the angle between S3 and S4 being taken 
as# — e gives the angle 29q between S2 and S4. Hence 
the first two clusters begin to show the spiral on the odd 
sites and the spiral on the even sites. Continuing this to 
the next few clusters shows that this allows a description 
of a system with two spirals, one on the odd, the other 
on the even sites, both with the same turn-angle 28o or 
wavelength. Also the energy of each cluster is the same, 
so one can consider just one cluster, say the first one 
above, in figuring out the relation between e and d. For 
the case of 28q = ir/2, spins order as following pattern 
(\t/ 1 — >\iy / < — \). Therefor, in the case (6^0), 
the following solutions can satisfy the minimum energy 
condition, 

[6,9') = (0 > 0),(7r,7r) J (7r,0) ) (0,7r) J 

[6,9') = (6>o,-0o + e), (Spiral- type), where 

168a 2 - (1 - 5) 2 
C ° S0 ° ~ ~Aa(l-Sni + SY 
sin 6*o = — >• 9q = nil — > e = 0. 

(14) 
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By substituting them into the Eq.©, the ground state 
energy of cluster in different sectors becomes 

ft-antifcrro = h c (lT, 7r) = S 2 ( - 1 + a) , 

Kud = h c (0, 7r) = S 2 {-6 - a), 
hdvm = h c (n,0) = S 2 (S - a), 
h _ u I a a , ,\ _ o2 / 



^spiral 



= h e (e ,-e + e) = s' 



1-6 2 
46 ■ 

J 2 — 1 
45 ' 



(15) 



These critical lines separate spiral phase from the uud 
and the duu phases. There is a discontinuity in the 
derivative ^gg-, and therefore a first-order phase tran- 
sition through the mentioned critical lines. In addi- 
tion, one should note that the ordering of the uud 
and duu phases, in principle is a type of the stripe- 
antiferromagnetic phas o 21 ' 22 . Therefore, the order pa- 
rameter for distinguishing these phases is defined as 



M: 



N 



E(- 1 )" +1 ( 5 2„-l+^in), ^ UUd 



71 = 1 



_ J2(-l) n+1 (S^ + Sl n+1 ), for duu.(16) 



We have also found antiferromagnetic phase that is stable 
for 

l ~ 62 t X n 
a < ; — for 5 > 0, 



4 

<5 2 -l 



a > 



for 5 < 0. 



(17) 



It is completely clear that in uud and duu phases, M* p 
takes the value 1. In FIG. @] (a), we have plotted M* t 
as a function of the frustration a for different values of 
the frustration 6 = 0.4, 0.6, 0.8. As it can seen from this 
figure, for a < a Cl the staggered magnetization is equal 
to 1.0, which shows that the ground state of the system 
is in the fully polarized antiferromagnetic phase. Finally, 
it is clear that there is no antiferromagnetic ordering in 
the region a > a c . In FIG. S] (b), we have the same 
picture for the stripe antiferromagnetic order function. 



1.0 



0.5 - 



+ (1 + cos 6 ){ (1 - 6) + V(l - S) 2 cos 2 O + 45}) . 

As it can be seen from the above equations, in respect 
to the case of S = 0, dimerization exchange removed the 
degeneracy between (|, |, t) an d (t, |> 4-) states. The (I, t 
, t) state, defines as a phase with opposite magnetization 
on odd (J(l — S)) bonds, but (t, til) state denotes by 
opposite magnetization on even ( J(l + 5)) bonds. Using 
the conditions in Eq. (Ti"4"|) allow us to find the stability of 
different phases. Doing some calculations, one can find 
two critical lines as 
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FIG. 4: (Color online). The staggered magnetization M„t and 
stripe-antiferromagnetic M 3p order parameters as function of 
the frustration a for different values of the frustration 5 = 
0.4,0.6,0.8. 



The zero value of M* p in the region a < ag is in com- 
plete agreement with the fully polarized antiferromag- 
netic and spiral phases in this region. By more increas- 
ing the dimerization and for a > ag, clearly be seen that 
the ground state of the system is in the uud. Also one 
predicts that the tripe antiferromagnetic as a function 
of the a displays a jump for certain parameters which is 
one of the most important indications of the first-order 
phase transition. We emphasize that a first-order phase 
transition occurs between the spiral and the duu phases, 
for negative values of the dimerization. 

The FIG. ((5]) shows the exact classical ground state 
phase diagram of the model in 5 — a plane. It should 
be mentioned that the same phase diagram can be also 
found in the ferromagnetic side which we do not depict. 
In the absence of dimerization, 5 = 0, there are antiferro- 
magnetic (ferromagnetic) and spiral phases which is sepa- 
rated by two critical point at a c = ±0.25 which the neg- 
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FIG. 5: (Color online). The classical ground state Phase 
diagram of the dimerized and frustrated chains in 5 — a plane. 
The antiferromagnetic is separated from spiral by critical line. 
At 5 = 0, the spiral phase is separated from antiferromagnetic 
phass by the critical point a = 0.25. 



We check the possible configurations which can mini- 
mize the above equations. By taking arbitrary ((/>, </>'), 
we have the solutions (tt, tt), (0, 0) that are related to 
the antiferromagnetic and ferromagnetic states respec- 
tively. We have also the solutions (tt, 0) — > (I, f, t) and 
(0, tt) —> (t, f, I) that are related to the uud and duu 
states respectively. Spiral phase also exists same as the 
coplanar case. The stability of different phases in the 
non-coplanar case is also checked and behaves same as 
the coplanar case. 



ative sign refers to ferromagnetic side. The second order 
phase transition occurs at these critical points. By turn- 
ing dimcrization the spiral phase persists to be stable in 
i^p 1 < a < region for 5 > and ^^-4 < a < ^ 1 
region for 5 < 0. The antiferromagnetic phase remains 

1 —A 2 

stable up to the critical line ° . The uud and duu 
phases propagate with different energy and separated by 
two first order critical lines, ( and S 4S 1 ), from spiral 
phase. 

In the following we are interested to implement LK 
method to non-coplanar antiferromagnetic case. Again 
with using cluster approach and without losing gener- 
ality, labeled the angles (9,<fi) ,(9',<fi') made by the end 
spins with the central spin (see FIG. [2]). Then the cluster 
energy is determine as 



h c (e,e') = 5 2 {( 



1-5 



cos 9 



1 + 5 



cos 9' 



cos 9 cos 9 



' + sin 9 sin 9' cos(4> - j . 



(18) 

Minimizing over the angles 9,9',4> , and 4>' gives the 
following equations 



dhc 
89 

dhc 
89' 

dhc 
d<j> 
dhc 



C-2 . 

1(1 -5) sin 6 + 2a [sin 9 cos 9' 

cos 9 sin 9' cos(</> — </>')]} , 
S 2 r 

-y|(l + <5) sin0 / + 2a[sin0 / cos (9 

cos 9' sin6*cos((/) — </>')] j, 
S 2 



| - 2a sin 9 sin 9' sin(0 - (j)') j , 
^- |2a sin 9 sin 0'sin(0 -</>')}■ (19) 



2 

S 2 



III. CONCLUSION 



To summarize, we have studied the classical ground 
state magnetic phase diagram of the dimerized and 
frustrated Heisenbcrg chain using LK cluster method. 
In coplanar case and in the absence of dimcriza- 
tion effect this approach could detect antiferromag- 
netic (ferromagnetic) and spiral phases. We have shown 
that turning the dimcrization yields to remove the degen- 
eracy between two uud and duu phases. We have argued 
that in the ground state phase diagram of the system 
there are first order transition lines. These lines separate 
spiral and uud or duu phases. On the other hand two sec- 
ond order phase transition points also exist, which sepa- 
rate antiferromagnetic(ferromagnetic) and spiral phases. 
By helping this approach we have calculated the spiral 
exact critical exponent e = 1/2. 



In order to generalize our treatment, we have consid- 
ered the non-coplanar case and checked its phases by LK 
cluster method. Our calculations revealed that in the 
non-coplanar case classical phase diagram consist of anti- 
ferromagnetic or ferromagnetic depend on nearest neigh- 
bor coupling, duu and uud phases which are still non- 
degenerate and spiral phase. Finally, one of the main 
achievement of our work which should be highlighted is 
that in the both coplanar and non-coplanar cases, the 
spiral state is stable in the classical phase diagram for 
5 ^ 0. 
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